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We analyze electron transport through relatively short coherent conductors in the presence 
of Coulomb interaction. We evaluate the current- voltage characteristics of such conductors 
taking into account the effect of an external environment. Within our model, at large con- 
ductances and low T the conductance is suppressed by a universal factor which depends only 
on the type of the conductor. We also argue that at T — the system "scatterer+shunt" 
can be either an insulator or a metal depending on whether its total resistance is larger or 
smaller than Rq — h/e 2 w 25.8 kf2. In a metallic phase the Coulomb gap is fully suppressed 
by quantum fluctuations. 



1 Introduction 

Discrete nature of the electron charge plays a crucial role in various phenomena in mesoscopic 
physics, causing, for instance, Coulomb blockade of electron tunneling in metallic junctions BE! 
and shot noise in mesoscopic conductors It is of particular importance to understand how an 
interplay of charge discreteness, Coulomb interaction and coherent scattering may affect electron 
transport in disordered conductors at sufficiently low temperatures. 

Recently we argued H that the interaction term in the current-voltage characteristics of a 
(relatively short) coherent conductor is controlled by the parameter 

P= EnT " (1 ~ Tn \ (1) 

which is already well known in the theory of shot noise @. Here T n are the transmissions of 
conducting modes. A similar conclusion was also reached in Ref. 6 in the limit of a single 



conducting mode, in which case the parameter 
In this paper we extend our previous results! 



]D reduces to (3 = 1 — T\. 

by considering the effect of an external environ- 



ment on quantum transport through a coherent scatterer in the presence of Coulomb interaction. 
2 Quasiclassical Langevin equation 

As in Ref. 5 we will consider electron transport through an arbitrary coherent scatterer between 
two big reservoirs. The scatterer conductance without interactions is defined by the standard 
Landauer formula 1/R= (2e 2 /h) J2 n T n . Phase and energy relaxation are only allowed in the 
reservoirs and not during scattering, i.e. the scatterer is shorter than both dephasing and 
inelastic relaxation lengths L„ and Lj n . Coulomb effects in the scatterer region are described by 
an effective capacitance C. We also assume that the scatterer is attached to an external voltage 
source V x via a linear impedance Zg(u). 

In the limit of sufficiently large energies and/or at large scatterer conductances g = Rq/R 3> 
1 it is convenient to describe the system by means of the quasiclassical Langevin equation 
approach. This equation can be derived from the real time path integral technique BSi, For a 
coherent conductor to be considered here this derivation was performed in Ref. 5. The effect of 



a linear external impedance can be incorporated in the same way as it was done for the case of 
tunnel junctions^!. Combining the results!'! one arrives at the following Langevin equation 
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Here <p(t)/e = V(t) is the fluctuating voltage across the conductor. Eq. (|2|) is sufficient provided 
eV and T are smaller than the typical inverse traversal time 1/rtrav (e.g. the Thouless energy 
in the case of diffusive conductors). Below we will mainly address the limit RC > r trav , however 
some of our results should remain valid in the opposite case as well. The terms in the right-hand 
side of (^) account for the current noise. They are defined by the correlators 
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If we set cp(t) = eV x t and define the total fluctuating current 6I(t) = £1 cos eV x t+£,2 sineV^t+ 
£3, we will immediately reproduce the standart result 0'!: 
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The I — V curve for a conductor can be obtained by averaging eq. (H|) over noise. We find 

IR = V-([£icoB(p + &Ban(p]). (6) 

The last term in eq. (||) describes the effect of Coulomb interaction. We note that this term 
depends only on the two stochastic variables £1 and £2. Since the correlation functions for both 
these variables (|3|) are proportional to the parameter /3, the magnitude of the whole interaction 
term in <M) should scale with the same parameter. Thus, the result (0) takes the form 
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where f(V,T) is the universal function which depends on R and Zgioj). This function was 
already evaluated in the case of tunnel junctionsB. Defining \/Z{u) = 1/R — wC + l/Zg(u) 
and proceeding perturbatively in ReZ, at T — > one finds (cf., e.g., Refs. 9 and 6) 
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In a special case of a linear Ohmic environment Z${u 
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where go = Rq/Rq = g + gs and gs = Rq/Rs- Eqs. fli~0| ) work well provided either go S> 1 
or max(T, eV) 3> Ec- In the limit go S> 1 and max(eV,T) ^> goEc exp(— go/2) one can set 
exp(-F(t)) ~ 1. Then eqs. (0), (|) yield the result 
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where ^(x) is the digamma function, w 

dl 



u + iv, u = goEc/n 2 T and v = eV/irT. At T 



from (11) we obtain 



R 



dV 



1 



go ln i 1 + (eVR Cy 



(12) 



while in the limit eV/Ec 2> max(l, g) we find RI = V — (3e/2C. At V — > from (|iT|) we get 
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(13) 



The above logarithmic dependencies on eV and T should also hold for RqC < T tra v, in which 
case in eqs. (12) and (|l3| ) with the logarithmic accuracy one can set goEc — > 1/Tt mu . 

At very low temperatures and voltages max(T, eV) < goEc exp(— go/2) the interaction cor- 
rection to the conductance of a coherent scatterer saturates at a universal value which does not 
depend on the interaction but only on the transmission distribution. This result follows imme- 
diately from Eqs. (§), (0)- Evaluating (0) at long times we find ~ (2/g )(ln(t/R C) +7) 
and performing the integral in (||), in the leading order in /?/<?o "C 1 we obtain 
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This formula successfully reproduces a complete Coulomb blockade of tunneling G — > in the 
limit /3 — ► 1 (tunnel junctions), demonstrates the absence of it for ballistic scatterers {(5 — > 0), 
and yields suppression of the Landauer conductance by the factor 2/3 for diffusive conductors. 



3 Comparison with experiments 

In order to compare our theoretical results with recent experiments we will use the data reported 
recently by two groups I1M3. 

Weber et al. I1U experimentally investigated the I — V curve of a short (L ~ 90 nm L v , Lj n ) 
diffusive conductor fabricated as a bridge between two big metallic reservoirs. Even though the 
conductance of this bridge was large g ~ 2000 the Coulomb blockade effect on the I — V curve 
was clearly visible becoming more pronounced with decreasing temperature, see Fig. 3 of Ref. 
11. It was observed that the experimental datalHI were well described by the formula 



G(V, T) = G(0,T ) + A\n(T/T ) +A$(eV/T), 



(15) 



where A was determined to be A ~ (0.4 0.7) x Rq 1 depending on the sample and $(eV/T) 
was found to be a universal function of eV/T which tends to zero at V — * 0. The logarithmic 
dependence of G(0, T) on temperature (|l5|) agrees with our eqs. (@), (pi). Both the function 
$>(eV/T) and the value A can easily be evaluated within our theory. Making use of eq. ( |i~T| ) in 
the experimentally relevant interval of temperatures and voltages we obtain A = 2(3 /Rq and 
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The function (|i~6| ) is plotted in Fig. 1 together with the experimental data0 for the sample No. 
1. One observes a very good agreement between theory and experiment except at high voltages 
where heating effects become dominant. From the same comparison for the sample No. 1 one 
finds P ~ 0.25 slightly smaller than the value (5 = 1/3 expected in a diffusive limit. 

In another experiment Krupenin et al. 113 studied the I — V curves of Cr resistors fabricated 
in the f orm of 2d strips and also clearly observed Coulomb blockade effects. The ~ 1/xm long 
samples 113 were not diffusive (their resistances were too large to be described by the Drude 
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Figure 1: The function (|l6| ) plotted together with the experimental data for the sample No 1. Different 
experimental curves correspond to different temperatures. 
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Figure 2: Eq. (0) plotted together with the experimental dataEJ for the sample f with R ~ 2.5 Kf2 per square. 



formula) and most likely had a granular structure. For our comparison we will adopt the model 
of a granular array with grains connected via short coherent scatterers with some average value 
of the parameter /3. For the sampled with the resistance 2.5 kf2 per square one can estimate the 



average conductance of one scatterer as g ~ 10. Again, a very good agreement between theory 
and experiment is observed, see Fig. 2. From the best fit one can extract the value ~ CL3 5. 

A relatively low value of the scatterer conductance g ~ 10 reached in the experiments!^ also 
allows to enter the regime max(T, eV) < gEc exp(— g/2), where the conductance was predicted 
to saturate at the level (|l4|). According to this prediction the ratio between the values of G 
measured in the limits of low and high voltages should approach 1 — j3. From the data 111 we 
extract f3 ~ 0.31. This value is close to one found from the fit in Fig. 2. Thus the datal!3 clearly 
support both our theoretical predictions ( |l2|) and (|i~4|). It is also interesting that both estimates 
give the values of f3 very close to (3 = 1/3 which one would expect for diffusive scatter ers. 

Finally, we notice that the results of recent conductance measurements 1121 on multi-wall 
carbon nanotubes, with typical lengths of order L« or shorter, are also consistent with the 
logarithmic dependence predicted in eq. ([12]). Further experimental work on the subject would 
be highly desirable. 



4 Instantons and metal-insulator phase transition 



Under which conditions does the result (14) remain valid? One of them was already formulated 
above: go ^> 1, i.e. either g or g$ should be much larger than one. In this case the conductance 
G saturates at the value ( |i~4| ) for max(eV, T) < goEc exp(— go/2). In order to establish another 
important limitation for the result (|l4]) let us recall that, even though our quasiclassical Langevin 
equation approach does account for nonlinear in (p effects, this approach may nevertheless become 
insufficient at very low energies even for g$ ^> 1 because it does not include instantonsi-3, i.e. 
nontrivial saddle points for the exact effective action. 

An important part of the instanton analysis was carried out by Nazarov £3 who derived 
the renormalized Coulomb energy Eq for a general coherent conductor within the exponential 
accuracy. The result can be written in the form Ec oc Ec exp(— ag). One can also go beyond 
the exponential accuracy and estimate the pre-exponent in the expressioiL for Ec- This was 
done in Ref. 16. Without going into details here, we only quote the result £3 
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This formula is valid for ag ^> 1, i.e. either at large conductances g 2> 1 or, if g ~ 1, for very 
small values R n implying o> 1. 

The quantity Ec (0) plays the role of an effective Coulomb gap in our problem. At T < Ec 
our Langevin equation analysis is insufficient and, hence, eq. (14) becomes inaccurate. In this 
regime the conductance decreases with T and the system is an insulator at T = 0. 

Let us now take into account the effect of an external linear Ohmic environment Z$ — Rs- 
In this case quantum fluctuations of the charge in the shunt resistor will affect the Coulomb gap 
Ec and - as we shall see - may even lead to its total^unpression provided Rs is sufficiently low. 

In order to proceed we will follow the analysis EJ>£9. Treating the charge q as a quantum 
variable □ and integrating out the phase tp one can map our problem onto that of a linearly 
damped quantum particle q in a periodic potential. This is a well-known problem0'i which can 
be treated, e.g., by means of the renormalization group (RG) technique. Successively reducing 
the high frequency cutoff uj c and integrating out charges with higher frequencies one arrives 
at the standard RG equations £3. After trivial manipulations these equations can be rewritten 
directly in terms of the combination ag = 5Z n ln R^ 1 . One obtains 



d(ag)/d(hxu} c ) = (1 - gz) (1 + l/ag), dgx/d(lnuj c ) = 0, 



(18) 

where gs = ggs/9o is the total dimensionless conductance. Eqs. ( |l8| ) are valid as long as ag ^> 1. 
One observes that for 5s < 1 the quantity ag decreases in the course of renormalization. Hence, 



in that case the Coulomb gap remains nonzero, the charge q is localized and the system is an 
insulator at T = 0. The effective Coulomb gap Ec can be defined as the energy scale at which 
the renormalized value ag becomes of order one. Then from ( |i~8| ) one finds 

i 

E c ~ Eciagexpi-ag)] 1 -^ . (19) 

On the other hand, for gs > 1 the combination ag always scales to larger values. In this case the 
Coulomb gap is fully suppressed by quantum fluctuations, Eq = 0, the charge q is delocalized 
and the conductance remains nonzero ([l4|) even at T = 0. This is a metallic phase. A quantum 
phase transition between the insulating and metallic phases occurs at g-£ = 1- 

Finally, let us briefly discuss possible implications of our results for recent experiments 
which strongly indicate the presence of a metal-insulator phase transition in various 2d disor- 
dered systems. One can consider a (sufficiently small) coherent scatterer with the dimensionless 
conductance g viewing all other scatterers in the system as an effective environment with the 
conductance gs- Assuming this environment to be Ohmic at sufficiently low frequencies, one 
immediately arrives at the conclusion about the presence of a quantum metal-insulator phase 
transition at gs=l. In 2d systems one has g ~ gs ~ go ~ <?£. Therefore in such systems this 
phase transition should be expected at conductances ~ 1/ Rq , exactly as it was observed in many 
experiments EJ. Local properties of the insulating and metallic phases are expected to be very 
different. In the insulating phase charges should be localized around inhomogeneities (puddles) 
due to Coulomb blockade, while in the metallic phase the Coulomb gap is suppressed and the 
charge distribution should be much more uniform. These expectations are fully consistent with 
experimental observations El. Thus, there might be a direct relation between the experimental 
resultsli3li3 and the old problem of a dissipative quantum phase transition til. 
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